RIORDAN MATRIX REPRESENTATIONS OF 
EULER'S CONSTANT 7 AND EULER'S NUMBER e 



EDRAY HERBER COINS AND ASAMOAH NKWANTA 
This paper is dedicated to David Harold Blackwell (April 24, 1919 - July 8, 2010). 

Abstract. We show that the Euler-Mascheroni constant 7 and Euler's 
number e can both be represented as a product of a Riordan matrix and 
certain row and column vectors. 



1. Introduction 
It was shown by Renter |2j that the Euler-Mascheroni constant 
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can be represented as a product of an infinite-dimensional row vector, the 
inverse of a lower triangular matrix, and an infinite-dimensional column 
vector: 
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Renter's proof uses induction, definite integrals, convergence of power se- 
ries, and Abel's Theorem. In this paper, we recast this statement using 
the language of Riordan matrices. We exhibit another proof as well as a 
generalization. Our main result is the 
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Theorem 1. Consider sequences {ao, a\, . . . , a n , . . . }, {bo, bi, . . . , b n , . . .} 
and {co, ci, . . . , c n , . . . } of complex numbers such that ao, bo, Co / 0; as well 
as an integer exponent d. Assume that 

(i) the power series a(x) = J2 n a n x n , b(x) = J2 n b n x n , c(x) = Y, n c n x n , 
and b(x) d are convergent in the interval \x\ < 1; and 

(ii) the following complex residue exists: 



Res 



z=0 



a(z)b(z~ 1 ) d c(z- 1 ) 



-L/ a(z)b(z- 1 ) d c(z- i )^. 

2 ~Kl J\ z \ = l Z 



Then the matrix product 



(a ai a 2 



62 



60 
61 



b 



b n 



bo 



Cl 
C2 



/ W 



is eguaZ to the above residue. 



The infinite-dimensional lower triangular matrix is an example of a Riordan 
matrix. Specifically, it is that Riordan matrix associated with the power 
series b(x) d . Renter's result follows by careful analysis of the power series 
(1.3) 

log(l-x) 1 1 1 2 1 

= = l + -x+-x 2 + ■■■ + — -x n + --- 

x 2 3 n + 1 



a(x) 
b(x)- 1 
c(x) 



x 



1 



1 



1 



log (1 - x) 
a(x) — 1 



~ 1 2 X 12 X 24 X "' LnXU 

1 1 1 2 1 n 

« + n X + 7 3 ' + ■ ■ ■ + ~ X™ + • • • 

2 3 4 n+2 



The coefficients L„ are sometimes called the "logarithmic numbers" or the 
"Gregory coefficients" ; these are basically the Bernoulli numbers of the sec- 
ond kind up to a choice of sign. (Renter employs the coefficients Ck = —L^.) 
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The idea of this paper is that we have the matrix product 



(1.4) 
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24 



V : 7 



: ••• : \ : / 

which is equivalent to the recursive identity 53m=o ^m/ (n — m) = 0, that is 
valid whenever n = 2, 3, 4, . . . . The matrix product, and hence the recursive 
identity, can be derived from properties of Riordan matrices. Renter's result 
follows from the identity Ylm=i L m /m = 7, which in turn follows from an 
identity involving a definite integral. 

As another consequence of our main result, we can also show that Euler's 
number 



(1.5) 



lim 

n— too 



1 

1 + - 

n 



2.7182818284. 



can be represented as a product of an infinite-dimensional row vector, a 
lower triangular matrix, and an infinite-dimensional column vector. 

Corollary 2. For any integers p, q, and d with pq > 1, the number 

dn 



pq 



pq 



— v e d = lim 

1 n^rco 



pq 



pq — 1 
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is equal to the matrix product 
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In the process of proving these generalizations, we present a representation 
theoretic view of Riordan matrices. That is, we consider the matrices as 
representations ir : G — > GL(V) of a certain group G - namely, the Riordan 
group - acting on an infinite-dimensional vector space V - namely, the 
collection of those formal power series h(x) in C[x] where h(0) = 0. 
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2. Introduction to Riordan Matrices 



We wish to list several key results in the theory of Riordan matrices. To 
do so, we recast this theory using techniques from representation theory 
very much in the spirit of Bacher [TJ. Our ultimate goal in this section is to 
explain how Riordan matrices are connected to a permutation representation 
7r : G — > GL(V) of a certain group G acting on an infinite dimensional 
vector space V. Some of the notation in the sequel will differ from standard 
notation such as that give by Shapiro et al. [4j and Sprugnoli [5], [6], but 
we will explain the connection. 

2.1. Group Actions. Before developing the representation theoretic view, 
we give the definition of a Riordan matrix and few related useful properties. 
Let A; be a field; it is customary to set k = C as the set of complex numbers, 
but, in practice, k = Q is the set of rational numbers. Set as the 

collection of formal power series in an indeterminate x; we will view this as 
a A;- vector space with countable basis {l, n , . . . }. For most of 

this article, we will not be concerned with regions of convergence for these 
series. 

There are three binary operations kfx} x k^xj —> k^xj which will be of 
importance to us, namely multiplication •, composition o, and addition +. 
Explicitly, if we write 

oo oo 



n=0 



(2.1) f( x ) = Y J fnX n and g( x )=Y J 9n 

n=0 

then we have the formal power series 

n 

^ ] fm 9n- 



X 



(2.2) 



(/•<?)(*) = E 

n=0 

oo 

[fog)(x) = J2 



if + 9) 
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oo 
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oo 
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fn + 9n 



There are three subsets of the vector space k{x\ which will be of interest to 
us in the sequel. 



Proposition 3. Define the subsets 



H 



K 



V 



fix) e klxj 
g[x) e klxj 



h(x) G k{xj 



g(0) = yet g'(0) ^ 



h(0) = 
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(i) H is a group under multiplication •, K is a group under composition 
o, and V is a group under addition +. In particular, V is a k-vector 
space with countable basis {x, x 2 , . . . , x n , . . . }. 

(ii) The map ip : K — > Aut(H) which sends g(x) G K to the automor- 
phism (fig : f(x) 1 y (f °g)(x) is a group homomorphism, where g(x) 
is the compositional inverse of g{x). In particular, G = H K is 
a group under the binary operation * : G x G — > G defined by 

(fi, #i)*(/2, 92) = (fi • <P 91 (f2), 9i og 2 ). 

(iii) The map * : G x V — > V defined by (/, g)*h = f • (ho~g) is a group 
action of G onV . 

We use g(x) to denote the compositional inverse g" 1 {x) so that we will 
not confuse this with the multiplicative inverse g(x)~ 1 . Later, we will show 
that G is isomorphic to the Riordan group R. Moreover, we will show that 
H, a normal subgroup of G, is isomorphic to the Appell subgroup of R. 
The motivation of this result is to use the action of G on V to write down 
a permutation representation tt : G — > GL(V), then use the canonical basis 



{ 



\ . . . } of V to list infinite-dimensional matrices. 



Proof We show (i) to fix some notation to be used in the sequel. Since 
(/•s)(0) = /(0) g(0) / for any /(x), g(x) G H, we see that • : HxH H 
is an associative binary operation. The identity is the constant power series 
e(x) = 1, and the inverse of /(x) is its reciprocal, seen to be a power series 
by expressing as a formal geometric series: 



fix) 



(2.3) 
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Jo 



CO 
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n=0 



Y,i-fn/h)x n 



n=0 

00 
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m=0 m+' 
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■+n m = 
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Since (/o <?)(()) = f(g(0)) = /(0) = and (/o^)'(O) = f'(g(0))g'(0) = 
f'iO)g'(0) / for any f(x), g(x) G i^T, we see that o : K x i^T ->■ K is an 
associative binary operation. The identity is the power series id(x) = x, and 
the inverse of g(x) is its compositional inverse g(x) = Y^ n ~g n x n having the 
implicitly defined coefficients 

9o = 0, 
1 

9i 



(2.4) 



91 



E 

m=0 



.9,; 



E 



9, 



5n 







for n = 2, 3, 



Since (/ + #)(0) = /(0) + 5(0) = for any /(x), 5 (x) G V, we see that 
+ : V x V — > V is an associative binary operation. The identity is the 
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constant power series o(x) = 0, and the inverse of h(x) is the negation 
—h(x), seen to be a power series with (— h)(0) = —h(0) = 0. 

Now we show (ii). Since (/ o g)(0) = f(g(0)) = /(0) + for any f(x) G 
H and g(x) G K, we see that <p : K — > Aut(iT) is well-defined. Given 
h(x) £ K we have <p g oip h = ip goh because for all f(x) £ H we have 



V 9 °^ [/0*0] = V?<? 



(foh)(o 



(2.5) 



= ^/ ° h o g ) (x) = [ fogoh)(x) 
= <Pgoh[f(x) 



Hence <p : K — > Aut(-ff) is indeed a group homomorphism. The semi-direct 
product G = H K consists of pairs (f{x), g(x)) with f(x) G i? and 
g{x) G if, where the binary operation * : G x G — > G is defined by 



(2.6) /!(*), gi(x))*(f 2 (x), g 2 (x)) = ( fi(x) f 2 (gl(x)), gifaix)) 



Finally, we show (iii). The map * : G x V — > V is defined as the formal 
identity 



f(x), g(x))*h(x) = f(x)h(g(x)). 



(2.7) 



Since [(/, g)*h](0) = /(0) h(g(Q)) = /(0)/i(0) = 0, we see that the map 
* : Gx V — > V is well-defined. As the identity element of G is (e(x), id(x)) = 
(l, x), we see that (e(x), id(x))*h(x) = h(x) so that it acts trivially on 
V. Given two elements (/i, g\), (f 2 , g 2 ) G G and /i(x) G F, we have the 
identity 



Slip), gi(x) * 



/2(ar), 52 (x) )*h(x) 



/i(x), 5i (x) * 



(2.f 



/2(ar) h(g 2 (x)) 

fi(x) f2(gi(x)) h(g^ogj(x)) 
fi(x) S2{gT{x)) h{g 1 og 2 {x)) 

Si(x) f 2 (gl{x)), g 1 (g 2 (x)))*h(x) 



Si(x), gi(x) * / 2 (z), g 2 (x 



*h(x). 
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Similarly, given two elements h\{x), h,2(x) G V and (/, g)& G, we have the 
identity 



f(x), g(x) * 



h\(x) + /l2(ic) 



(2.9) =f(x) h!(g(x)) + fo(g(x)) 

= g(x)j*hi(x)+ (f(x), g(x) j*h 2 (x). 

Hence * : G x V — > V is indeed a group action. □ 
2.2. Riordan Matrices. Recall that the set 



(2.10) 



V = \ h(x) G kfxj 



h(Q) = 



is a k- vector space with countable basis [x, x 2 , . . . , x n , . . . }. Since the 
semi-direct product G = H ys^ K acts on V, we have a "permutation" 
representation it : G — > GL(V). Explicitly, this representation is defined on 
the basis elements of V via the formal identity 



(2.11) 



f(x),g(x))*x m = f(x)[g(x)Y 

00 

= ^ ^ ln,m X 



for m = 1, 2, 3, ... . 



n=l 



(Recall that g(x) is the compositional inverse of g(x).) The matrix with 
respect to the basis {x, x 2 , . . . , x n , . . . } is given by the lower triangular 



matrix 



(2.12) 



7T f{x), g(x) 



h,i h,2 h,s 



\ 



l>n,l ln,2 l>n,3 ' ' ' l>n,n 

V i : ; ... : 



■7 



Recall that g(0) = yet /(0), </(0) 7^ 0. The following result explains the 
main multiplicative property of these matrices. 

Theorem 4. Continue notation as above. 

(i) 7r : G — > GL(V) is a group homomorphism. That is, 



n[fi(x), gi(x)\ 7T (/ 2 (x), S2(x)l =tt /i(x)/ 2 (5i(x)), 51(52(2;)) 
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(ii) For a generating function t(x) = to + t\ x + t 2 x 2 + • • • with to 7^ 0, 



) j — ( ^n,p ip-m J 
' Vp=i / n , m >l 



Such matrices 7r(/, 5) are called the Riordan matrices associated to the 
pair (/, g) . The collection R of Riordan matrices is a group which is isomor- 
phic to G = H y\tp K; this is the Riordan group. The collection of matrices 
7r(/, id) is a group which is isomorphic to H; this normal subgroup is the 
Appell subgroup of R. 

Proof. We show (i) . In the proof of Proposition [3l we found that for each 
h{x) € V we have the following formal identity involving power series as 
elements of k fxj : 



fi(x), gi(x) * 



h(x), g 2 (x) )*h(x) 



(2.13) 



fi(x), gi(x) )*[f 2 (x), g 2 {x 



th(x) 



h(x) f 2 (gi(x)), gi(g 2 {x)) ) *h(x) 



In particular, this holds for the basis elements h(x) = x n , so the result 
follows. 

Now we show (ii). For a generating function t(x) = to + 1\ x + 1 2 x 2 + • • • , 
we have the product 

00 



n—m i 



(2.14) (t(x), id(x)\ *x m = t(x) x m = f2 t 

so that matrices in the Appell subgroup are in the form 

h to 



x ; 



(2.15) 



7r( t(x), id(x) 



t 2 h to 

tn-l t n - 2 t n -3 

\ ■ ■ : 

This gives the matrix product 

n(f(x), g(x)\ n\t(x), x^j = (l njP 

(2.16) 



to 



7 



n,p>l 



'p—m 



p,m>l 



m \ 
^ ^ ln,p tp—m J 
P =l / n,m>l 
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so the result follows. 



□ 



2.3. Examples. Let k = Q. Using elementary Calculus, we find the power 
series expansions 



(2.17) 



ln(l-x) , 1 1 2 1 3 1 4 1 5 

i = l + - x + -x 2 + -x d + -x 4 + -x b + 

x 2 3 4 5 6 



x 



1 1 1 9 1 3 19 4 3 5 
1 X X X X X + 

2 12 24 720 160 



ln(l-x) 

which are valid on whenever \x\ < 1. Hence the formal power series 
(2.18) /(x ) = i + I x + I x 2 + l x 3 + ... + _^_ x „ + ... 

is an element of H, and has multiplicative inverse 
1 



(2.19) 



fix) 

We have the product 



1 1 1 2 1 3 19 4 3 5 
1 X X X X X + 

2 12 24 720 160 



(2.20) (f(x) t id(x)) *x™ = /(x) x™ = f; ra _^ + 1 

^ ' n=l 



X 



which yields the matrix 



(2.21) 



7T(/, id) = 



/l 

1 

2 

1 

3 



I _J_ ... 1 

n n— 1 n— 2 



Similarly, we have the product 
1 



(2.22) 



-, id(x) *x 



= x m — - x m+l 3- x m+2 3- x m+3 — ^"^ x m+4 + 

2 12 24 720 
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Since we may use Theorem U] to conclude that 7r(/, id) 1 = 7r(l//, id), we 

find the identity 

(2.23) 

r 1 



i 

2 

1 

3 



1 

n-1 



V 



7 



1 




1 

2 


1 


1 


1 


12 


2 


1 


1 


24 


12 


19 


1 


720 


24 



1 

2 

J_ 

' 12 



\ 



These matrices are elements of the Appell subgroup of R. 



2.4. Relation with Standard Notation. Standard references for Riordan 
matrices are Shapiro et al. [I] and Sprugnoli [5], [6]. The notation 7r(/, g) 
employed above is not the typical one, so we explain the connection. Con- 
sider sequences {G , Gi, G 2 , . . . , G n , ...} and {Fx, F 2 , F 3 , . . . , F n , . . .} 
of complex numbers k = C, where Go, F\ ^ 0. Upon associating generating 

functions G(x) = G +Gi x+G 2 x 2 -\ and F(x) = F t x+F 2 x 2 +F 3 x 3 -\ 

to these sequences respectively, the standard notation for a Riordan matrix 
is that infinite-dimensional matrix given by 



(2.24) 



L= [G(x), F(x)] =ir(G(x), F(x)) = (l 



m 

n,m>l 



in terms of the compositional inverse F(x) of F(x). Indeed, the entry l n , r 
in the nth row and mth column satisfies the relation 



3& 



(2.25) G(x)[F(x)] m = ^2l n , m x n for m = 1, 2, 3, ... 

n=l 

as formal power series in C[a?]. Equivalently, a Riordan matrix L can be 
defined by a pair (G(x), F(x)) of generating functions. 

Corollary 5 (Fundamental Theorem of the Riordan Group: [3], [4], [6]). 

Continue notation as above. 

(i) The product of Riordan matrices is again a Riordan matrix. Explic- 
itly, their product satisfies the relation 



[Gx(x), Fx(x)] [G 2 (x), F 2 (x)] 



Gx(x)G 2 (Fx(x)), F 2 (Fx(x)) 
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(ii) For a generating function T(x) = Tq + T\ x + T 2 x 2 + • • • with To ^ 0, 
we have the product 



[G(x), F(x)] [T(x), 



x 



VP=1 



p—m 



n,m>l 



Proof. Statement (i) is shown in [H Eq. 5] and O Proof of Thm. 2.1], but 
we give an alternate proof. Upon denoting fi(x) = Gi{x) and gi{x) = Fi(x) 
for i = 1 and 2, we find the matrix product 

[GiOc), [G 2 (x), F 2 (s)] 



(2.26) 



7T 



ffi(^)) 7r^/2(ai), 52(2;) 
fi(x) f2{gT(x)), 9\{92{x)) 

G l (x)G 2 {F l (x)), F 2 {F 1 (x)) 

which follows directly from Theorem [H Statement (ii) is also shown in [3], 
but it follows directly from Theorem 0] as well. □ 

3. Proof of Renter's Result and Generalizations 
3.1. Main Result. We now prove the following: 

Theorem[T], Consider sequences {ao, a\, . . . , a n , ■ ■ ■}, {bo, bi, . . . , b n , . . . } 

and {cq, ci, . . . , c n , . . . } of complex numbers such that ao, bo, cq 7^ 0; as well 
as an integer exponent d. Assume that 

(i) the power series a(x) = Yj n a n x n , b[x) = £] n b n x n , c(x) = ^ n c n x n , 
and b{x) d are convergent in the interval \x\ < 1; and 

(ii) the following complex residue exists: 



Res 



2=0 



a(z) b(z 



- l ) d c(z- v ) 



-L/ a(z)b(z- 1 ) d c(z- 1 )^. 
2vrz /| 2 | = i z 



Then the matrix product 



fb 

h 
b 2 



(ao «i a 2 



b 
h 



\ fco\ 

Cl 

C 2 



b n b n -\ b n -2 



bo 



■7 W 



is equal to the above residue. 
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Proof. With the three power series a(x) = ^2 n a n x n , b(x) = ^ n b n x n , and 
c(x) = c n x n convergent in the interval \x\ < 1, consider the power series 



(3.1) 



f(x) = b(x) d c(x) = ^2f n x n where |x| < 1. 



n=0 



As elements of the Appell subgroup of R, we invoke Theorem [4] to see 
that we have the matrix product 7r(/(x), x) = 7r(6(x), x) d 7r(c(x), x). In 
particular, the first column is given by 



(3.2) 



h 
h 

in 



/b 
h b 
b 2 h 



bo 



b n b n -\ 6 n _2 



b 



Cl 
C-2 



Hence the matrix product 
(3.3) 



a ai a 2 



(bo 
bi b 
b 2 h 



bo 



b n b n -\ b n -2 



Cl 
C2 



is equal to the sum ^ n a n f n . We wish to evaluate this sum using complex 
analysis. 

By assumption, the power series a(x), b(x), and c{x) are convergent in the 
interval \x\ < 1. Hence, for each fixed real number r satisfying < r < 1, 
the functions a(z) and f(z) are uniformly convergent inside a closed disk 
| z | < r. Hence we can interchange summation and integration to find the 



RIORDAN MATRIX REPRESENTATIONS OF 7 AND e 



13 



integral around the boundary to be equal to 



2ni 



— 6 a(z)b(z*) d c(z*) — 



dz 



\z\=r 



1 



dz 



— 6 a(*) /(**)- 



(3.4) 



27ri J\z\=r 

00 00 

y ; y ] a ni fn,2 r 

ni=0 ri2=Q 
00 

„2n 



+ri2 



2vr 



2tt 



3 i(ni-n 2 )fl ^ 



n=0 



Here 2* is the complex conjugate of z. As r — > 1, the integral exists so by 
Cauchy's Residue Theorem it must be equal to 

'a(z)b(z- 1 ) d c(z- 1 y 



Res 



2 = 



(3.5) 



lim 

r— >1 



lim 

r-+l 



2 

1 

27^ 



z|=l z 



^ ] fn 



„2n 



n=0 



^ ^ a n fn- 



n=0 



The Theorem follows upon equating this with equation (|3.3 



□ 

3.2. Applications. We explain how to use Theorem Q] in order to express 
Euler's number e = 2.7182818284 ... in terms of Riordan matrices. 

Corollary [2], For any integers p, q, and d with pq > 1, the number 

/ t \ dn 
P<1 lim ^-r(l + ± 



pq-1 

is equal to the matrix product 



pt 



1 



pn 



(I 



/1\ 
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Proof. The coefficients of the matrices above correspond to the three power 
series 



(3.6) 



a(x) 
b(x) 
c{x) 



1 — x/p 



1 

1 - x/q 



Ex 

n=0 F 



£ 

n=0 

oo 

£ 

n=0 



where \x\ < 1. 



For a complex number z with \z\ < 1, we have the identity 



(3.7) 



a{z)b{z- 1 ) d c{z~ l ) 


1 








1 


z 


1 — z/p 








1 - z- 1 /q_ 



£ 



E 



d™ 2 



,ni— 112— n'i=n+l 



p Ul n 2 ! q ns 



The residue corresponds to the coefficient of the z 1 term, so we consider 
the terms where n = — 1: 



Res 



z=0 



a(z)b(z' l ) d c(z- 1 y 



E 



d n2 



(3.8) 



ni=ri2+nz 
oo 



p n ^ n 2 ! q n3 



E 

,n 2 =0 

pq-1 



n 2 ! \p 



The Corollary follows now from Theorem [TJ 



□ 



Renter's result is also an application of Theorem [TJ 
Corollary 6 (|2j. The Euler-Mascheroni constant 



7 = lim 

n— too 



n ^ 



Inn 



Vm=l 



0.5772156649 . 



RIORDAN MATRIX REPRESENTATIONS OF 7 AND e 
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is equal to the matrix product 



(1 



/I 

1 

2 

1 

3 



1 1 

2 3 



r 1 / n 



1 1 1 

n n— 1 n— 2 



1 

n+1 



: '■■/ V : / 



Proof. The coefficients of the matrices above correspond to the three power 
series 



a(x) 

(3.9) b(x) 
c(x) 



log (1 - x) 

X 

log (1 - x) 

X 

a(x) — 1 



n=0 



ra=0 



ra + 2 



where Ixl < 1. 



We will choose the exponent d = — 1. We will express the reciprocal as the 
power series 

19 



= -1 + - x H x z + — + — x 4 H x° + 



log(l-x) " ' 2 ' 12 '24 ' 720 ' 160 
(3.10) 

n=0 

which is also convergent in the interval \x\ < 1. (Recall that the coeffi- 
cients L„ are sometimes called the "logarithmic numbers" or the "Gregory 
coefficients".) For a complex number z with \z\ < 1, we have the identity 



(3.11) 



q(z)6(z- 1 ) Q 'c(z- 1 ) 

z 

log(l - z) 



+ 



log(l - z)" 



= E^t+ E 



z 
00 

E 



log(l-z-i) 



n=0 



n + m + 1 



z . 
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The residue corresponds to the coefficient of the z 1 term, so we consider 
the terms where n = — 1: 



Res, 



a(z)b(z- 1 ) e c(z- L ) 



(3.12) 



OO j 



m=l 
1 



in 



o 



,m=l 



„m— 1 



1 1 

X log (1 - 



dx 



= 7- 

The Corollary follows now from Theorem [TJ □ 

We conclude with stating that Theorem [1] can also be used to show Rior- 
dan matrix representations for In 2 and 7r 2 /6. Finding matrix representations 
of other constants, like \/2, ir, and the Golden Ratio </>, are of interest. 
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